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Systematic calculations of fission barriers allowing for triaxial deformation are performed for even- 
even superheavy nuclei with charge number Z = 112 — 120 using three classes of covariant density 
functional models. The softness of nuclei in the triaxial plane leads to an emergence of several 
competing fission pathes in the region of the inner fission barrier in some of these nuclei. The outer 
fission barriers are considerably affected by triaxiality and octupole deformation. General trends of 
the evolution of the inner and the outer fission barrier heights are discussed as a function of the 
particle numbers. 

PACS numbers: 21.60.Jz, 24.75.+i, 27.90. +b 



o 

=5 
S3 



> 
00 

m 

in 
o 

(N 



I. INTRODUCTION 

The on-going search for new superheavy elements is 
motivated by the attempts to provide theoretical and ex- 
perimental explanations for two open questions in nu- 
clear structure. The first question is related to the limits 
of the existence of atomic nuclei at large values of proton 
number, while the second one to the location of the is- 
land of stability of superheavy nuclei and the next magic 
numbers (if any) beyond Z = 82 and N = 126. Heavy 
and superheavy nuclei decay by spontaneous fission and 
therefore their stability is defined essentially by the size 
and the shape of the fission barriers. Thus, the inner fis- 
sion barrier is considered a fundamental characteristic of 
these nuclei which is important in resolving the questions 
mentioned above. 

A systematic investigation of the properties of fission 
barriers is the best way to address the stability of nu- 
clei against spontaneous fission since it eliminates the 
arbitrariness of conclusions with respect to the choice 
of a specific nucleus. Different models with triaxiality 
included have been used for extensive theoretical inves- 
tigations of the properties of inner fission barriers in the 
actinides. It has been found that the heights of inner fis- 
sion barriers are reduced when triaxial shapes are taken 
into account However, this reduction has a strong 

dependence on the particle number in the proton and 
neutron subsystems as well as on the applied model. 
These investigations were performed in the following 
frameworks: microscopic+macroscopic (MM) methods 
(3411, the extended Thomas- Fermi plus Strutinsky inte- 
gral (ETFSI) method [l(j, and non-relativistic energy 
density functionals (EDF) based on zero-range Skyrme 
@, EMU and finite-range Gogny @, MM forces, and 
recently in covariant density functional theory (CDFT) 

These methods have also been used for the study of 
fission barriers in superheavy nuclei (see Table IIVI below 



for a review). However, only few of them take into ac- 
count triaxiality. For example, within the framework of 
covariant density functional theory, it has been consid- 
ered only in the stud y o f fission barriers in a single nu- 
cleus ( 264 Hs) in Ref. [llj. Thus, in order to fill this gap 
in our knowledge a systematic study of fission barriers 
of even-even Z = 112 — 120 superheavy nuclei with tri- 
axiality included is performed in the current manuscript. 
We use the same method that has been successfully cm- 
ployed in Ref. [l| for the study of inner fission barriers in 
actinides, where an average deviation from experiment 
of 0.76 MeV has been found. When extrapolating to 
superheavy nuclei we want to understand the impact of 
triaxiality on fission barriers and how they evolve with 
the change of the particle numbers. An additional goal 
is to see how the fission barriers depend on the choice of 
the specific CDFT model. While the systematic investi- 
gation of inner fission barriers of actinides in Ref. [l[ has 
been performed with the nonlinear meson-nucleon cou- 
pling model represented by the NL3* parametrization of 
the relativistic mean field (RMF) Lagrangian, in the cur- 
rent study we use in addition to the nonlinear coupling 
model also density-dependent meson-nucleon and point 
coupling models. 

Density functional theories are extremely useful for the 
microscopic description of quantum mechanical many- 
body systems. They have been applied with great success 
for many years in Coulombic systems [IH, [2(| , where they 
are in principle exact and where the functional can be de- 
rived without any phcnomenological adjustments directly 
from the Coulomb interaction. In nuclear physics with 
spin and isospin degrees of freedom, with strong nucleon- 
nucleon and three-body forces the situation is much more 
complicated. However, covariant density functionals ex- 
ploit basic properties of QCD at low energies, in partic- 
ular symmetries and the separation of scales [2l|. They 
provide a consistent treatment of the spin degrees of free- 
dom, they include the complicated interplay between the 
large Lorentz scalar and vector self-energies induced on 
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the QCD level by the in-medium changes of the scalar 
and vector quark condensates [22|. In addition, these 
functionals include nuclear magnetism [23j], i.e. a con- 
sistent description of currents and time-odd mean fields 
important for odd- mass nuclei [24|, the excitations with 
unsaturated spins, magnetic moments [25j and nuclear 
rotations [2H 127} . Because of Lorentz invariance no new 
adjustable parameters are required for the time-odd parts 
of the mean fields. Of course, at present, all attempts to 
derive these functionals directly from the bare forces [28l — 
l3lj | do not reach the required accuracy. However, in re- 
cent years modern and very successful covariant density 
functionals have been derived [321-135} which are based on 
density dependent vertices and one additional parameter 
characterizing the range of the force. They provide an 
excellent description of ground states and excited states 
all over the periodic table [3(1 [37} with a high predictive 
power. Modern versions of these forces derive the density 
dependence of the vertices from state-of-the-art ab-initio 
calculations and use only the remaining few parameters 
for a fine tuning of experimental masses in finite spheri- 
cal [sH and deformed |Hj] nuclei. 

The manuscript is organized as follows. The theoret- 
ical framework and the details of the numerical calcula- 
tions are discussed in Sec. [TTJ The results of the investi- 
gations of the fission barriers, the role of triaxiality and 
the comparison with the actinidc region are presented in 
Sec. IIIII Finally, Sec. IIVI summarizes the results of our 
work. 



II. THEORETICAL FRAMEWORK AND THE 
DETAILS OF NUMERICAL CALCULATIONS 



Three classes of covariant density functional mod- 
els are used throughout this manuscript: the nonlin- 
ear meson- nucleon coupling model (NL), the density- 
dependent meson-exchange model (DD-ME) and a 
density-dependent point coupling model (DD-PC). The 
main differences between them lay in the treatment of the 
range of the interaction, the mesons and density depen- 
dence. The interaction in the first two classes has a finite 
range, while the third class uses zero-range interaction 
with one additional gradient term in the scalar-isoscalar 
channel. The mesons are absent in the density-dependent 
point coupling model. The density dependence is explicit 
in the last two models, while it shows up via the non- 
linearity in the cr-meson in the nonlinear meson-nucleon 
coupling model. Each of these classes is represented in 
the current manuscript by a set of parameters that is 
considered to be state-of-the-art. 



is a standard Lagrangian density [4C 



In the meson-exchange models [33|, [38|, |39| , the nucleus 
is described as a system of Dirac nucleons interacting 
via the exchange of mesons with finite masses leading to 
finite-range interactions. The starting point of covariant 
density functional theory (CDFT) for these two models 



C = t/j [j(id - g u ui - g p p f - eA) - m- g a a\ ip 



-n pu w + x -ml^ (1) 



which contains nucleons described by the Dirac spinors 
ip with the mass m and several effective mesons char- 
acterized by the quantum numbers of spin, parity, and 
isospin. They create effective fields in a Dirac equation, 
which corresponds to the Kohn-Sham equation [l9| in 
the non-relativistic case. The Lagrangian ([T]) contains as 
parameters the meson masses m a , m u , and m p and the 
coupling constants 170-, g u , and g p . e is the charge of the 
protons and it vanishes for neutrons. This linear model 
has first been introduced by Walecka [4l|, [42j . 

To treat the density dependence in this model Boguta 
and Bodmer [43| introduced a density dependence via a 
non-linear meson coupling replacing the term ^m 2 a 2 in 
Eq. Q by 



U(a) 



1 



f 



1 



m 2 a 2 + -g20-° + -.g-iv" ■ 



(2) 



The nonlinear meson-nucleon coupling is represented by 
the parameter set NL3* [H[ (see Table [I]), which is a 
modern version of the widely used parameter set NL3 
pi} . Apart from the fixed values for the masses m, m w 
and m p , there are six phenomcnological parameters m a , 
ga, 9u, 9 P , 92, and g 3 . 

The density-dependent meson-nucleon coupling model 
has an explicit density dependence for the meson-nucleon 
vertices. There are no nonlinear terms in the a meson, 
i.e. <72 = .93 = 0. The meson-nucleon vertices are defined 
as: 



9i{p) = 9i(p8at)fi(x) for i = a,u),p 
where the density dependence is given by 

1 + h(x + di) 2 



fi{x) = at 



l + dix + diY 



for a and U) and by 



f p (x) = cxp(-a p (:r - 1)). 



(3) 



(4) 



(5) 



for the p meson, x is defined as the ratio between the 
baryonic density p at a specific location and the baryonic 
density at saturation p sat in symmetric nuclear matter. 
The eight parameters in Eq. ^ are not independent, 
but constrained as follows: fi(l) = 1, f a (1) = 
and /j (0) = 0. These constraints reduce the number of 
independent parameters for density dependence to three. 
This model is represented in the present investigations 
by the parameter set DD-ME2 33( given in Table HI 
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TABLE I: The parameters of the NL3* and DD-ME2 pa- 
rameterizations of the Lagrangian. Note that g a — gcr(psat), 
Puj = g^(p S at) and g p = ff p (p sa t) in the case of the DD-ME2 
parametrization. 



Parameter 


NL3* 


DD-ME2 


m 


939 


939 


m a 


502.5742 


550.1238 


m u 


782.600 


783.000 


m p 


763.000 


763.000 


9<T 


10.0944 


10.5396 




12.8065 


13.0189 


g P 


4.5748 


3.6836 


92 


-10.8093 




93 


-30.1486 




(la 




1.3881 


K 




1.0943 


Ca 




1.7057 


d a 




0.4421 






1.3892 


b„ 




0.9240 


Cu, 




1.4620 


(L 




0.4775 


a p 




0.5647 



TABLE II: The parameters of the DD-PC1 parameterization 
in the RMF Lagrangian 



Parameter 


DD-PC1 


m 


939 


a CT 


-10.04616 


b a 


-9.15042 


C<, 


-6.42729 


d a 


1.37235 


a-u 


5.91946 


b u 


8.86370 


d u 


0.65835 


b P 


1.83595 


dp 


0.64025 



TABLE III: The G", G%, G\ and G v 2 parameters [in MeV] for 
different parameterizations of the RMF Lagrangian. 

Force G" G£ G{ G p 2 

NL3* 10.7 -10.4 7.40 18.9 

DD-PC1 10.5 -7.38 7.5 19.2 

DD-ME2 11.0 -10.3 7.9 17.0 



The Lagrangian for the density-dependent point cou- 
pling model (34 . [45( 1 is given by 

C = ip (17 • d — to) ip 
1 



a s (p) (V>V>) {ipi>) - ^avip) (^7 M V>) {i>1ti<f) 
arv(p) (tpf^ip) (ipTj^tp) 



2 Ss (d v W>) ( W) ~ efa ■ A {1 2 T3 V - (6) 



It contains the free-nucleon Lagrangian, the point cou- 
pling interaction terms, and the coupling of the proton 
to the electromagnetic field. The derivative terms in Eq. 
© account for the leading effects of finite-range inter- 
action which are important in nuclei. In analogy with 
meson-exchange models, this model contains isoscalar- 
scalar, isoscalar-vector and isovector-vector interactions. 
In the present work it is represented by the DD-PC1 
parametrization 34| given in Table [TTJ 

The triaxial relativistic mean field (RMF)+BCS ap- 
proach [46[ is used here for the description of fission bar- 
riers. This approach has been very successfully applied 
to a systematic description of the fission barriers in the 
actinides The RMF-equations are solved iteratively 
and at each iteration the BCS occupation probabilities 
v\ are determined. These quantities are used in the cal- 
culation of densities, energies and new fields for the next 
iteration. We use a monopole pairing force [47]] with the 
strength parameters G T for neutrons (r = n) and protons 
(r = p). 

For each type of particle we start with a pairing 
strength parameter G and solve at each iteration the gap 



equation [47 



l - = Y — 



(7) 



with Ek = \J (£fc — A) 2 + A 2 , where Ek are the eigenval- 
ues of the Dirac equation and the chemical potential A 
is determined by the average particle number. Then the 
occupation probabilities 



1 - 



A 



and the gap parameters 



A = G ^2 u k v k 

k>0 



(8) 



(9) 



arc determined in a self-consistent way. The pairing en- 
ergy is defined as 



k>0 



(10) 



The sum over k in Eqs. 0, © and flTU|) run over all 
states in the pairing window Eh < E cuto tf. In Ref. (48j 
empirical pairing gap parameters 



4.8 

A emp = — — MeV 



a omp _ 
P Z l/3 



MeV (11) 



have been determined by a fit to experimental data on 
neutron and proton gaps in the normal deformed mini- 
mum. 
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P 2 cos(y+30) 



FIG. 1: (Color online) Potential energy surface of the 240 Pu 
nucleus with the NL3* parametrization of the RMF La- 
grangian. The energy difference between two neighboring 
cquipotential lines is equal to 0.5 MeV. The pink dashed line 
with solid circles shows the lowest-energy solution as a func- 
tion of 02 . Small red stars show the positions of deformation 
points at which the numerical results were obtained. The fig- 
ure is based on the results of Ref. p}. Further details are 
given in the text. 

These empirical gap parameters form the basis for the 
definition of the strength parameters G T in the current 
manuscript. The following expressions [5Cj| 

A • G n = — C?2 — 7— MeV (12) 

A-G p = G\ + G p 2 ^-j^ MeV (13) 

arc used in the calculations. First, using empirical 
gap parameters of Eq. (fTTjl . the values G n (Z,N) and 
G P (Z, N) are obtained for the ground states of all even- 
even nuclei in the Z = 112 - 126, N - Z = 48 - 62 
region. Then, the parameters G™, G%, G\ and G\ are 
defined by the least square fit to the set of the G n (Z, N) 
and G P (Z, N). Their values depend on the parameter set 
of the Lagrangian and they are given in Table Hill In this 
way we have strength parameters for the effective pairing 
interaction depending in a smooth way on the neutron 
and proton numbers and, because of the changing level 
density, the gap parameters derived from those values 
show fluctuations as a function of the particle numbers. 
Note that similarly to Ref. [![, the cutoff energy for the 
pairing window is set to E cuto ff=120 MeV. 

The calculations are performed by successive diagonal- 
izations using the Broyden method [58| and the method 
of quadratic constraints [47|. We have also implemented 
the augmented Lagrangian method [59| for constraints 
in our computer codes, but we did not find any clear ad- 
vantage of this method over the method of quadratic con- 
straints for the type of potential energy surfaces (PES) 



we are dealing with. Starting from the three multipolc 
operators 

Q 20 = 2z 2 -x 2 -y 2 (14) 
Q22 = x 2 -y 2 (15) 
Q 3 o = z(2z 2 - 3.t 2 - 3y 2 ) (16) 

we calculate in the following investigations three types 
of potential energy surfaces, using three types of con- 
straints. 

First, we restrict ourselves to axially symmetric config- 
urations with reflection symmetry, abbreviated through- 
out the paper by (A). Here we use the computer code 
DIZ [fjoj based on an expansion of the Dirac spinors in 
terms of harmonic oscillator wave functions with cylin- 
drical symmetry and we minimize 

(H) + C 20 ((Q 20 ) - q 20 ) 2 (17) 

where (H) is the total energy, (Q20} denotes the expecta- 
tion values of the mass quadrupole operators, q 2 $ is the 
constrained value of the multipole moment, and G20 the 
corresponding stiffness constant |47j . 

In the next step abbreviated throughout the paper by 
(T) we use the triaxial computer code DIC [23[ with the 
D2-symmetry based on an expansion of the Dirac spinors 
in terms of a cartesian oscillator basis and by imposing 
constraints on axial and triaxial mass quadrupole mo- 
ments we minimize 

(H)+ J2 C 2 ^{Q 2p )-q 2 rf. (18) 

^1=0,2 

In addition, starting from the axial reflection symmet- 
ric computer code DIZ of Ref. [6(|, we have developed 
an axial reflection asymmetric [octupolc deformed] code 
DOZ. It is used for a study of the impact of axial octupolc 
deformation on the inner and outer fission barriers. The 
calculations in this axially symmetric octupole code, ab- 
breviated throughout this paper by (O) vary the function 

(H)+ Y, C Mi o«Q Mi o)-?m,o) 3 . (19) 

M=2,3 

It was checked that the numerical results obtained for 
axially symmetric solutions in the A, T and O calcula- 
tions differ by no more than 50 keV for all deformations 
of interest. 

The truncation of the basis is performed in all these 
calculations in such a way that all states belonging to 
the shells up to Np = 20 fermionic shells and Nb = 20 
bosonic shells are taken into account. This truncation 
scheme has been tested and used in the actinides in 
Ref. [l|. Although the calculations in such a trunca- 
tion scheme provide sufficient numerical accuracy, they 
are also very computationally demanding (see Ref. [l[ for 
details) . This is the reason why we treat the pairing chan- 
nel in the present triaxial RMF calculations in the BCS 
approximation despite the fact that the triaxial cranked 
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FIG. 2: (Color online) The difference between experimental and calculated heights of inner fission barriers as a function of 
neutron number TV. The results of the calculations are compared to estimated fission barrier heights given in the RIPL-2 
database |4iJ, which is used for this purpose in the absolute majority of theoretical studies on fission barriers in actinides. The 
results of the calculations within microscopic+macroscopic method ('MM(Dobrowolski)' @] and 'MM (Moller)' Q), covariant 
density functional theory ('CDFT' pi) and density functional theory based on the finite range Gogny force ('Gogny DFT' 
[3) are shown. Thick dashed lines are used to show the average trend of the deviations between theory and experiment as 
a function of neutron number. The average deviation per barrier S [in MeV] is defined as 5 = £f =1 \B)(th) - B)(exp)\/N, 
where N is the number of the barriers with known experimental heights, and Bf(th) (Bj(exp)) are calculated (experimental) 
heights of the barriers. Long-dashed lines represent the trend of the deviations between theory and experiment as a function 
of neutron number. They are obtained via linear regression based on a least square fit. 
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FIG. 3: (Color online) The same as in Fig. [2] but as a function of proton number Z. 



Rclativistic Hartree+Bogoliubov (RHB) approach with 
the finite range Gogny forces in the pairing channel was 
developed in the late nineties EH • RMF+BCS cal- 
culations are less time-consuming and more stable (espe- 
cially, in the saddle point region) than the RHB calcula- 
tions. 

Fig. [T] shows an example of the distribution of defor- 
mation points on a potential energy surface (PES) at 
which the numerical results were obtained using these 



constraints in a triaxial (T) calculation. 

All mean-field calculations for fission barriers with tri- 
axiality included are computationally demanding. This 
is the reason why simple pairing interactions (either se- 
niority pairing with constant strength G or the zero-range 
(5-force) and the BCS framework are used in the majority 
of the calculations (see Table IIV|) . The overview given in 
this table focuses mainly on the investigations including 
triaxial deformation and on the most recent and com- 
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prehensive studies within the specific theoretical frame- 
works. In particular, the studies of fission barriers in 
superheavy elements (SHE) are included irrespective of 
whether the calculations include triaxial deformation or 
not. This table shows that the majority of studies neglect 
particle number projection and use different prescriptions 
for the size of the pairing window. The impact of these 
prescriptions on the height of the inner fission barrier has 
recently been investigated in Ref. [57| ■ 

The inclusion of triaxiality has improved the accuracy 
of the description of the inner fission barriers in actinides 
in all state-of-the-art models 0, @, 0, El- Figs. M and 
[3] show the differences between experimental and calcu- 
lated heights of inner fission barriers obtained in different 
theoretical models as a function of neutron and proton 
numbers, respectively. Note that this comparison covers 
only results of systematic triaxial calculations which si- 
multaneously include even-even Th, U, Pu, Cm and Cf 
nuclei. To our knowledge, no such calculations have been 
published with DFT based on Skyrme forces. As a result, 
these figures cover all existing systematic triaxial studies 
of inner fission barriers in actinides. 

The 8- values displayed on the panels of Figs. [5] and 
[3] show the average deviation from experiment for the 
calculated heights of inner fission barriers. One can see 
that they are of the same magnitude in the different ap- 
proaches and minor differences between the approaches 
in the 8- values are not important considering the con- 
siderable uncertainties in the extraction of inner fission 
barrier heights from experimental data as seen, for ex- 
ample, in the differences of the compilations of Refs. [HI 
and [651 ]. 

However, the similarity of the average trends of these 
deviations (shown by thick dashed lines in Figs. [5] and [3]) 
as a function of neutron and proton numbers is more im- 
portant considering the differences in underlying mean 
fields and in the treatment of pairing correlations. At 
present, it is difficult to find a clear explanation for these 
trends. Although differences in the treatment of pairing 
correlations (BCS with monopolc pairing and of differ- 
ent pairing windows in the CDFT [J and MM [|, 0] 
calculations versus the Hartree-Fock-Bogoliubov frame- 
work based on the D1S force in Gogny DFT Q) can 
contribute to deviations between theory and experiment 
[57j , it is quite unlikely that they are responsible for the 
observed trends of the deviations. 



III. RESULTS AND DISCUSSION 

The nucleus 292 120 is predicted to be a spherical dou- 
bly magic nucleus in CDFT [ll], [6(| . Its potential energy 
surface in the j3 — 7 plane is shown in Fig. @] It is inter- 
esting to compare it with the PES of the nucleus 240 Pu 
shown in Fig. [1] These two PES's are representative 
examples of typical PES's in actinides and superheavy 
nuclei. The gross structure of these two PES's is defined 
by the fact that the total energy is generally increasing 



when moving away from the 7 = 0° axis; so it looks like 
a canyon. However, there are local structures inside the 
canyon which define the differences between the two mass 
regions with respect to the impact of triaxiality on the 
inner and outer fission barriers. 

In 240 Pu, a large hill is located at the axial shape fa ~ 
0.5 inside a canyon. As a consequence, the fission path 
from the normal deformed minimum initially proceeds 
along the 7 = 0° axis, then bypasses the axial ,02=0.5 
hill via a path with 7 ~ 10°, and then proceeds along 
the bottom of the canyon on an axially symmetric path 
again. As a result of this bypass, the inner fission barrier 
heights of the actinides are lowered by 1 — 4 MeV due 
to triaxiality pj. However, the calculated outer fission 
barriers of the actinides are not affected by triaxiality 

Q. 

The properties of the PES of the nucleus 292 120, defin- 
ing the fission pathes, are completely opposite to the case 
of 240 Pu since there are two triaxial and one axial hills 
inside the PES canyon of 292 120. Two triaxial hills are lo- 
cated at moderate deformations (fa ~ 0.35, 7 ~ ±30°), 
while the axial hill is supcrdeformed (fa ~ 0.75). The 
fission path (shown by red dashed line in Fig. 2]) starts 
at a spherical shape, then proceeds between two triaxial 
hills (fa ~ 0.35, 7 ~30°) and bypasses the axial hill at 
fa -0.75 via a 7 - 7° path. The 7-softness of the PES, 
which exists between the two triaxial hills, has only a 
minor effect on the height of the inner fission barrier; the 
triaxial solution is lower than the axial one by 100-200 
keV at fa = 0.2 — 0.3 deformations (see Fig. [5]). However, 
this figure shows that along this fission path the height 
of inner fission barrier is not affected by triaxiality. 

The second fission path shown by a solid red line starts 
at spherical shape, and proceeds along the axially sym- 
metric 7 = 60° axis, via a saddle point at (fa — 0.47, 
7 ~ 21°) and then along the first fission path after second 
minimum. The existence of the valley between the walls 
of the canyon and triaxial hills at (fa — 0.35, 7 ~ ±30°) 
is a necessary condition for the existence of this path. 
This type of fission path exists in a number of nuclei, so 
for convenience we will call it as triaxial 'Tr-B' fission 
path. The unusual physical feature of this second path 
is the fact that initially the nucleus has to be squeezed 
along the axis of symmetry, thus creating an oblate nu- 
cleus with quadrupole deformation fa — —0.39. This is 
contrary to the usual picture of fission where the prolate 
or near-prolate nucleus is stretched out along the path of 
increasing quadrupole deformation. Note that the saddle 
of this fission path is only by 0.2 MeV lower than the sad- 
dle of the first fission path (see Table [Vj) and in addition 
it is much longer. Therefore, considering a dynamical 
calculation taking into account the action integral along 
the entire fission path [67| this path will probably not 
contribute much to the fission probability. 

Contrary to the actinides, the triaxiality has a consid- 
erable impact on the shape and the height of outer fission 
barrier which is lowered by ~ 3 MeV in 292 120 nucleus 
(see Fig. O. Note, that in this nucleus, the lowering of 



7 



outer fission barrier due to octupole deformation is sub- 
stantially smaller than the one due to triaxiality. 

It is also evident that the landscape of the PES and 
the existence of saddle points and valleys depends on the 
proton and neutron numbers. This is clearly seen in Fig. 
[Hand in Table fVl In the Z = 120 isotopes, the increase of 
neutron number up to TV = 184 and beyond leads to the 
emergence of an axial hill at fa ~ 0.2 — 25; this is clearly 
visible in the nucleus 304 120 in Fig. Its existence in 
high-Z and neutron-rich nuclei leads to the shift of the 
first fission path in the deformation space. In lighter sys- 
tems this path proceeds via the axially symmetric saddle 
as, for instance, in the nucleus 292 120 in Fig.|4] For heav- 
ier systems it starts at spherical shape 1 , proceeds along 
the axially symmetric 7 = 60° axis up to fa ~ 0.25 and 
then via a saddle point at (fa ~ 0.32,7 ~ 26°), located 
between a triaxial (fa ~ 0.35, 7 ~ 28 ) 2 and an axial 
fa ~ 0.2 hill, to the axially symmetric supcrdeformcd 
minimum. This type of fission path exists in a number of 
nuclei. We will label it as triaxial 'Tr-A' for convenience. 
In addition, there is a second triaxial path 'Tr-B' via a 
saddle at fa ~ 0.5,7 ~ 22°, which is similar to the one 
present in the nucleus 292 120. 

The structure of the PES of the nucleus 300 116 is simi- 
lar to the one of the nucleus 304 120 in Fig.SJ The lowering 
of neutron number below TV = 184 in the nuclei under 
study leads to the situation that the fission path via the 
axially symmetric inner saddle (called as 'Ax' below) be- 
comes energetically favored. This is clearly seen in the 
PES of the 286 116, 278 112 and 290 112 nuclei in Fig.Hand 
in Table El 

It is clear that the landscape of the PES of superheavy 
nuclei in the region of the inner fission barrier is more 
complicated than in the case of the actinides. The en- 
ergetically favored fission path in the majority of these 
nuclei proceeds through an axially symmetric saddle (see 
Fig. Hand Table Only in a few nuclei (304,306,308 120) 
300,302,304 118 and 300 116), triaxial fission path 'Tr-A' is 
lower in energy than axial 'Ax' fission path. The triax- 
ial fission path 'Tr-B' exists in the majority of the nuclei 
apart from neutron-poor Z = 114 nuclei and Z = 112 
nuclei (see Fig. [4] and Table lv|). The energy of the sad- 
dle point of this fission path is lower than the one of the 
first fission path (either 'Ax' or 'Tr-A') in a number of 
nuclei. However, the results of calculations within the 
MM method [13] suggest that although triaxiality lowers 
the static fission barriers, it plays a minor role in spon- 



1 As discussed later, some of the nuclei have superdeformed ground 
states. However, the calculated outer fission barriers are rather 
small (~ 2 MeV relatively to the superdeformed minimum) which 
suggests that these states are extremely unstable against fission. 
Thus, we start the discussion from sphcrical/wcakly(normal)- 
deformed minima which have better survival probability against 
fission. 

2 This hill is clearly visible when the energy difference between 
two neighboring equipotential lines is set to 0.1 MeV, see Fig. 1121 
below. 



taneous fission of superheavy nuclei with Z < 120. This 
is because a fission path via an oblate shape and triaxial 
saddles is substantially longer as compared to the axially 
symmetric path which leads to a significant reduction of 
the penetration probability. However, this complexity 
of the PES and the presence of two fission pathes for 
the inner fission barriers calls for finding the dynamical 
path along which the fission process takes place in CDFT. 
Work in this direction is in progress. 

Figs. [5] and [7] show deformation energy curves for sev- 
eral isotopes with charge numbers Z = 112, 114 and 116 
for three classes of CDFT models. Experimental esti- 
mates of the inner fission barrier heights were obtained 
for these nuclei in Rcf. (68|. The potential energy struc- 
ture of these nuclei is similar to the one seen in 292 120 
(Fig. 2]); the presence of two triaxial hills at moderate 
deformations (fa ~ 0.35, 7 ~ ±30°) defines two possible 
fission pathes: one between the hills (cither 'Ax' or 'Ax- 
Tr' pathes) and another between the hill and the walls 
of the PES canyon ('Tr-B' path). The 7-softness of the 
PES, which exists between the two triaxial hills, defines 
whether the fission path passing between these two hills 
is fully axially symmetric ('Ax' path) or has some de- 
gree of triaxiality ('Ax-Tr' path). The 'Ax-Tr' path is 
generally similar to the 'Ax' path but differs from it by 
the fact that moderate triaxiality (7 < 10°) appears at 
the saddle and along the shoulder of inner fission barrier. 
Although the heights of the discussed hills in the region 
of inner fission barrier depends on the parametrization, 
the observed features of the potential energy surfaces and 
fission pathes are rather independent of the parametriza- 
tion (Fig. |8j). Tables [V] and [VH] show that the saddle 
point along this path has triaxiality only in the case of 
the two nuclei 292 > 294 H4. However, even in these cases 
the 'Ax' saddle is only 50-150 keV higher than the 'Tr-A' 
saddle. Thus, similarly to the 292 120 nucleus triaxiality 
only marginally affects the inner fission barriers. Note 
that the ground states are somewhat deformed in these 
nuclei (see Figs. |5] and [§h). 

For these nuclei triaxiality has a considerable impact 
on the shape and height of the outer fission barriers; 
the decrease of the heights of the outer fission barri- 
ers due to triaxiality is typically in the range of 1.5-2.0 
MeV and it depends on the particle number and on the 
RMF parametrization. Note that the outer fission bar- 
riers are affected also by octupole deformation. In the 
Z = 112 nuclei, the triaxial saddle is lower in energy 
than the octupole saddle. The situation is reversed in 
the Z= 114, 116 nuclei. 

Among the different classes of CDFT models, the DD- 
ME2 parametrization always gives the highest values for 
the inner and outer fission barrier heights. They are (on 
average) by 1 and 1.5 MeV higher than the ones obtained 
in the NL3* and the DD-PC1 parametrizations. The 
heights and the shapes of the inner fission barriers are 
very similar for the NL3* and DD-PC1 parametrizations. 
The outer fission barriers also come close to each other 
in these two parametrizations in axial [A], triaxial [T] 
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and octupolc [O] (for Z = 112 nuclei) calculations. How- 
ever, with increasing of Z from 114 to 116 the difference 
between the results of octupole calculations increases; 
the difference between the energies of the octupole sad- 
dles obtained in the NL3* and DD-PC1 parametrizations 
reaches 2 MeV in the Z = 116 nuclei (see Fig. |6]). 

There is only one experimental work [68| where esti- 
mates on the heights of inner fission barriers in super- 
heavy nuclei with Z = 112, 114, and 116 have been ob- 
tained. Unfortunately, experimentally the fission barri- 
ers are accessible only indirectly and a model-dependent 
analysis is used to obtain these quantities, which causes 
an ambiguity in the comparison with theoretical results. 
Even in the actinide region where the fission barrier 
heights were extracted from a number of independent 
experiments with high statistics (see, for example, Ref. 
0), a typical uncertainty in the experimental values, as 
suggested by the differences among various compilations, 
is of the order of ±0.5 MeV (H)]. These uncertainties 
are expected to be higher in superheavy nuclei since the 
estimates of Ref. [68[ are based on experimental data rep- 
resented by low statistics and on a method which differs 
from the methods used in the analysis of fission barrier 
heights in actinides. In addition, there is no indepen- 
dent confirmation of the inner fission barrier height es- 
timates of Ref. HH . The interpretation of experimental 
data based on cross sections in terms of fission barrier 
heights becomes even more complicated when the fission 
path has a double hump structure, which according to 
many calculations may be the case in superheavy nuclei. 
The widening of the barrier due to the second hump (or 
its remnant) would require the lowering of the inner fis- 
sion barrier height; this possibility has not been taken 
into account in the analysis of Ref. j6^|. Based on this 
discussion it is clear that the level of confidence of fission 
barrier height estimates for superheavy nuclei is signifi- 
cantly lower than the one for the actinides. 

According to Ref. [68|, the estimated lower limits for 
fission barrier heights in even-even Z = 112, 114 and 116 
nuclei shown in Fig. [5] and [7] are 5.5, 6.7 and 6.4 MeV, 
respectively. Our results for the heights of inner fission 
barrier in these nuclei along the 'Ax' and 'Ax-Tr' fission 
pathes are always smaller than the experimental data by 
1 - 3 MeV. The saddle along the 'Tr-B' fission path is 
somewhat lower in energy (by 0.4-0.7 MeV) than the sad- 
dle along the 'Ax-Tr' path in ™ 2 ^ 4 llQ nuclei (see Table 
IVIip . However, as discussed above (based on the results 
of Ref. [63 ) it is unlikely that the fission predominantly 
proceeds along the 'Tr-B' path considering especially (i) 
that the gain in saddle point energy as compared with 
the 'Ax-Tr' path is small and (ii) the 'Tr-B' path is sub- 
stantially longer as compared with the 'Ax-Tr' path. 

Considering the discussion above, it is not clear at this 
moment how serious the discrepancy between calculations 
and experimental estimates is. However, the results of 
the calculations suggest one possible way to increase the 
heights of inner fission barriers. Potential energy surfaces 
in the ground state region of these nuclei are extremely 



soft (see Figs. [6] and [7]). For such nuclei, the correlations 
beyond mean field taken into account, for example, by 
the generator coordinate method can lower the energy of 
the ground state by a few MeV without affecting much 
the barrier top, thus effectively increasing the height of 
inner fission barrier. 

Systematic calculations of inner fission barriers have 
been performed for even-even Z = 112 — 120 nuclei with 
N — Z = 52 — 68 using the NL3* parametrization. This 
is the same parametrization which has been used in sys- 
tematic calculations of fission barriers in actinides [l| and 
which provides 0.76 MeV average deviation from exper- 
iment for the heights of inner fission barriers in these 
nuclei. For a given proton number Z a sequence of nine 
even-even nuclei is selected such that the middle nucleus 
of this sequence roughly corresponds to an experimen- 
tally observed nucleus (or a nucleus which is expected to 
be observed in near future). This crudely outlines the 
region which may be experimentally studied within the 
next one or two decades. The results of these calculations 
are summarized in Table fVl 

Based on the results of Ref. [63], it is reasonable to 
neglect the path 'Tr-B' since it is substantially longer as 
compared with other pathes which leads to significant 
reduction of the penetration probability. In Fig. [3Jd we 
consider the evolution of the heights of the inner fission 
barriers as a function of the neutron number TV. The 
Z = 112, 114 and 116 isotope chains show a generally 
increasing trend for the barrier heights with increasing 
neutron number for N > 172. Fig. [5^, suggests that the 
origin of this trend can be traced back to the deformation 
of the ground state. For small neutron numbers these 
nuclei are deformed in the ground state. However, they 
gradually become spherical when approaching N = 184 
because there is a spherical shell gap at this neutron num- 
ber (see, for example, Fig. 28 in Ref. 66]). The negative 
shell correction energy at the ground state is larger in 
absolute value in the vicinity of the N = 184 spheri- 
cal shell gap than at lower neutron numbers, where the 
ground state is deformed and characterized by a larger 
level density in the vicinity of the Fermi level. The level 
density (and, as a consequence, the shell correction en- 
ergy) at the saddle point of the inner fission barrier does 
not change so drastically as the one at the ground state. 
As a result, the heights of inner fission barriers, which 
are defined as the energy differences between the binding 
energies of the ground state and saddle point, show the 
observed features. Note that in these nuclei the energies 
of the 'Ax' and 'Ax-Tr' saddles differ by at most 200 keV, 
the difference between the energies of the 'Ax' and 'Tr-A' 
saddles does not exceed 400 keV. 

The nuclei with Z = 118 and 120 have, with the excep- 
tion of the nucleus 288 118, spherical minima due to the 
presence of the Z = 120 spherical shell gap. The nucleus 
292 120 has the highest value for the fission barrier among 
all nuclei under investigation. This is connected with its 
doubly magic nature in CDFT. In the isotope chain with 
Z = 120, moving away from the N = 172 shell closure, 
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shell effects connected with spherical shape become less 
pronounced and this leads to a decrease of the inner fis- 
sion barrier, because the barrier height in these nuclei is 
defined with respect to the spherical ground state. Apart 
from the lightest two and the heaviest Z = 118 isotopes, 
the fission barrier heights of the Z = 118 chain are nearly 
constant as a function of neutron number and they are 
close to 4.5 MeV. Note that in these nuclei the energies of 
the 'Ax' and 'Ax-Tr' saddles differ only by 100-400 keV. 
However, the difference between the energies of the 'Ax' 
and 'Tr-A' saddles can reach 2 MeV (as in the case of the 
nucleus 308 120, see Table E}. 

It is important to mention that the valley between the 
axial f3 2 ~ 0.2 and triaxial (3 2 ~ 0.40,7 ~ 25° hills lead- 
ing to a 'Tr-A' saddle is rather shallow. Its depth with 
respect to the tip of the triaxial hill varies between 100 
and 300 kcV. The latter value is obtained, for example, 
in 302 120 nucleus, see Fig.[T2l 

It is interesting to compare the current results with the 
ones obtained in other models. The results of Skyrme 
DFT calculations of Ref. [l]| for the N = 184 isotones 
show that the impact of triaxiality on the inner fission 
barrier is small for Z = 1 12, but increases with increas- 
ing Z (see Fig. 4 in Ref. [l3[ ). The lowering of the in- 
ner fission barrier due to triaxiality is around 2 MeV for 
Z = 120 and exceeds 3 MeV for Z = 126. In the ETFSI 
model calculations fioj the inner fission barriers are low- 
ered due to triaxiality in the nuclei (Z = 112, N = 182) 
and (Z = 114, N = 184) by 0.5 and 1.1 MeV, respec- 
tively. In the macroscopic+microscopic calculations of 
Ref. [|[ the largest reduction of the inner barrier height 
due to triaxiality is about 2 MeV and it appears in the 
region around Z 122, N « 180 (see Fig. 2 in Ref. [|). 

As discussed in Ref. [l[ on the example of actinide nu- 
clei, the reduction of the inner fission barrier height due 
to triaxiality is caused by the level densities in the vicin- 
ity of the Fermi level which are lower at triaxial shape as 
compared with axial one. The different location of the 
"magic" shell gaps in superheavy nuclei in the macro- 
scopic+microscopic model (at Z = 114, N = 184), in 
Skyrme DFT (predominantly at Z = 126, N = 184) and 
in CDFT (at Z = 120, N = 172) results in different 
deformed single-particle structures at the deformations 
typical for the saddles of the inner fission barrier. This 
is one of the sources of the differences in the predictions 
of different models. 

We have to keep in mind, however, that we use in this 
investigation a seniority zero pairing force with a fixed 
cutoff energy of E'cutoff = 120 MeV. It has been shown in 
Ref. [57| that pairing correlations play an important role 
for the calculation of the fission barriers and that even 
if the strength of the pairing force is adjusted to exper- 
imental gap parameters at the ground state, the range 
of the pairing force has an influence on the height of the 
barrier. This means for zero range forces and for the se- 
niority zero forces that the barriers depend on the cutoff 
energy. It also has been shown in Ref. [57| that the pa- 
rameter set DD-ME2 in connection with the finite range 



Gogny force D1S in the pairing channel produces in ax- 
ial symmetric calculations in most of the actinide nuclei 
inner barriers which are too high as compared with ex- 
periment (see Fig. 9 of Ref. [57|) and that it produces 
for the superheavy elements with Z = 112, 114, and 116, 
where experimental estimates are available [68j in axially 
symmetric calculations barriers which are in rather good 
agreement with those values (see Fig. 10 of Ref. [571]). 
Of course, so far, there exist no relativistic triaxial cal- 
culations with the finite range Gogn y fo rce in the pairing 
channel and full Gogny calculations [lj] are hard to com- 
pare because of the different spin-orbit force used in this 
model. We have, however to keep in mind, that triaxial- 
ity reduces the barriers in the actinides, but not in the 
superheavy elements with Z = 112, 114, and 116. There- 
fore, at the moment, we can only conclude that for a final 
comparison with the experimental data we have to wait 
for full RHB calculations or at least RMF+BCS calcu- 
lations with the finite range Gogny force in the pairing 
channel. Investigations in this direction are in progress. 

It is well known that reflection-asymmetric (octupolc 
deformed) shapes become important at the second fis- 
sion barrier and beyond (see Refs. 0, 0] and references 
therein). Our calculations indicate that triaxiality can 
play a similarly important role at the outer fission barri- 
ers. Usually, the triaxiality of these barriers is not men- 
tioned in the publications. To our knowledge, it is only 
discussed in Ref. @ that non-axial degrees of freedom 

play an important role in the description of outer fis- 
sion barriers of actinides. Fig. [TU] compares the energies 
of outer barrier saddle points as obtained in axial re- 
flection symmetric (A), triaxial (T) and axial reflection 
asymmetric [octupole deformed] (O) calculations. One 
can see that the inclusion of triaxiality or octupole de- 
formation always lowers the outer fission barrier. The un- 
derlying shell structure clearly defines which of the saddle 
points (triaxial or octupole deformed) is lower in energy. 
For example, the lowest saddle point is obtained in triax- 
ial calculations in proton-rich nuclei (N < 174). On the 
contrary, the lowest saddle point is obtained in octupole 
deformed calculations in neutron-rich nuclei (N > 174). 
Note that the decrease of the saddle point energy by oc- 
tupole deformation or triaxiality reaches 3 MeV in some 
nuclei. Thus, one can conclude that due to the struc- 
ture of the PES in the fission path valley, we observe in 
the superheavy region an opposite situation as compared 
to actinide nuclei where the triaxiality has no impact on 
the outer fission barriers. Our results also suggest that 
in some superheavy nuclei [mostly in the nuclei where 
the deformation energy curves of the A and O calcula- 
tions are similar in energy] the combination of two de- 
formations (triaxiality and odd-multipolc deformations) 
may be important in the definition of the fission path for 
/?2 > 0.5. The CDFT calculations with both deforma- 
tions included are at present not yet possible, but require 
further investigations. 

Fig. [TT] shows that some nuclei are superdeformed in 
the ground state. A summary for such nuclei is given in 
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Table IVTlTl Whether these states are stable or metastable 
should be defined by the height and the width of the 
outer fission barrier. The current calculations show that 
in many nuclei this barrier is appreciable in the axial cal- 
culations of type A. However, the inclusion of triaxial or 
octupole deformation decreases this barrier substantially 
so it is around 2 MeV in the majority of the nuclei. This 
low barrier would translate into a high penetration prob- 
ability for spontaneous fission, such that most likely these 
superdeformed states are metastable. Calculations of 
the spontaneous fission half-lives from spherical/weakly- 
deformed and superdeformed minima are needed in order 
to decide which of these minima is more stable against 
fission. Table |Vni] shows that the saddle point obtained 
in octupole deformed calculations (O) is the lowest in en- 
ergy in most of the cases. Only in proton-rich 294 120 and 
290,292-Qg nuclei, the saddle point of triaxial calculations 
(T) is the lowest. 

A superdeformed minimum exists also in the doubly 
magic nucleus 292 120 at a low excitation energy of ap- 
proximately 0.6 MeV (sec Fig. [5]). This is definitely 
not an artifact of the model under consideration, since 
similar minimum exists also in axial rclativistic Hartrcc- 
Bogoliubov calculations with the DD-ME2 parametriza- 
tion with Gogny PIS forces in the pairing channel (see 
Fig. 8 in Ref. |57[). Its calculated excitation energy de- 
pends on the actual strength of pairing and varies be- 
tween and 2 MeV. 

It is interesting to compare our results for the struc- 
ture of the outer fission barriers with those obtained in 
the axial RMF+BCS model (with and without reflec- 
tion symmetry) of Ref. [52[ which employed the NL3 and 
NL-Z2 parametrizations. Similar to our calculations, the 
superdeformed minima exist in the calculations of Ref. 
[52j (see Fig. 5 in this reference) without octupole defor- 
mation. However, in the calculations of Ref. [52[ with 
the NL-Z2 parametrization, the heights of outer fission 
barriers with respect to the SD minimum are lower by 
approximately 2 MeV as compared to our calculations. 
As a consequence, the inclusion of octupole deformation 
completely eliminates the outer fission barriers leading to 
results contradicting ours. On the other hand, the results 
of the calculations of Ref. [52| (see Fig. 6 in this reference) 
for a few selected nuclei based on the NL3 parametriza- 
tion are very similar to ours, namely, the ground state is 
superdeformed and the outer fission barrier has a height 
of approximately 2.5 MeV. This result is not surprising 
considering that the NL3* parametrization is very simi- 
lar to NL3 38]. However, the results of the calculations 
for a few selected Z = 112, 114 and 116 nuclei shown in 
Figs.[6]and[7]rcvcal that the outer fission barriers survive 
in the presence of octupole and triaxial deformation not 
only in the NL3* parametrization but also in the DD- 
ME2 and DD-PC1 parametrizations. On average, the 
height of outer fission barrier in these nuclei is around 2 
MeV. On the contrary, the calculations with the NL-Z2 
parametrization in Ref. [52j show that octupole defor- 
mation kills the outer fission barriers in these nuclei (see 



Fig. 5 in this reference). 

IV. CONCLUSIONS 

We have carried out first systematic investigations 
of fission barriers in even-even superheavy nuclei with 
Z = 112 — 120 within covariant density functional theory 
including triaxial shapes with D2-symmetry and octupole 
shapes with axial symmetry. Three different classes of 
models with the state-of-the-art parameterizations NL3*, 
DD-ME2 and DD-PC1 were used in the calculations. 
Pairing correlations are taken into account in the BCS 
approximation using seniority pairing forces adjusted to 
empirical values of the gap parameters. The following 
conclusions have been obtained: 

• The low-Z and low- AT nuclei in this region are char- 
acterized by axially symmetric inner fission bar- 
riers. The increase of the particle numbers leads 
to a softening of the potential energy surfaces in 
the triaxial plane. As a result, several competing 
fission pathes in the region of inner fission bar- 
rier emerge in some of the nuclei. Their impor- 
tance in spontaneous fission can be defined only 
by taking into account the fission dynamics more 
seriously. However, the results of the calculations 
within the macroscopic+microscopic method (67j 
suggests that the 'Tr-B' path (and even maybe 'Tr- 
A' one) may not be so important since they are 
substantially longer as compared with the axially 
symmetric path because it leads to a significant re- 
duction of the penetration probability. 

• Triaxiality lowers the outer fission barriers by 1.5-3 
MeV in reflection symmetric calculations. In many 
nuclei the lowering due to triaxiality is even more 
important than the one due to octupole deforma- 
tion, which known to be important for the outer fis- 
sion barriers and beyond from previous calculations 
(see Ref. [52j and references therein). The under- 
lying shell structure clearly defines that the triax- 
ial [octupole] saddle is lower in energy for proton- 
rich nuclei with N < 174 [neutron-rich nuclei with 
N > 174]. 

• On average, inner and outer fission barriers ob- 
tained for the NL3* and DD-PC1 parametriza- 
tions are similar. On the contrary, the DD-ME2 
parametrization produces barriers which are by 1- 
1.5 MeV higher than the ones obtained with NL3* 
and DD-PC1. 

• The superdeformed minimum is the lowest in en- 
ergy in some of these nuclei. In the present cal- 
culations, the outer fission barriers with respect to 
these minima are about 2 MeV high. Both minima 
and barriers are present in all three classes of the 
CDFT models. It has to be investigated, however, 
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if these superdcformcd minima are stable with re- 
spects to more general deformations not taken into 
account so far. 

The comparison of our results with those of non- 
rclativistic models clearly shows that CDFT predictions 
for the heights of inner fission barriers in the superheavy 
region with a seniority zero pairing force still remain on 
the lower end among nuclear structure models used so 
far. Only axially symmetric calculations with the finite 
range Gogny force D1S in the pairing channel can repro- 
duce the estimates of inner fission barrier heights of Rcf . 
[68j . Considering the uncertainties of these estimates, 
the investigation of other experimental observables that 
strongly depend on fission barrier heights, especially fis- 



sion half-lives and /3- or electron capture delayed fission, 
is needed. The work in this direction is in progress. 
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TABLE IV: The definition of pairing in the studies of inner fission barriers within the last decade. The first column 
shows the first author and the reference. The pairing model (BCS or HFB) and type of pairing (G - seniority 
pairing with fixed strength G, 8 - zero range 5-force, Gogny - finite range Gogny force) are shown in column 2. 
Column 3 shows either the region of the nuclear chart or the nucleus in which the fitting of the parameters of the 
pairing force has been performed in the case of constant G and/or 5-pairing. Column 4 shows whether particle 
number projection (PNP) by means of the Lipkin-Nogami method has been used in the calculations. Column 5 
indicates whether the systematic calculations, covering all even-even actinide nuclei with measured inner fission 
barriers, have been performed (Yes) or not (No). In the case of restricted calculations, the number of nuclei, 
for which the calculations have been performed, is shown. Column 6 is similar to column 5 but for SHE. The 
calculations are considered to be systematic if they cover a significant range of proton and neutron numbers. "T" 
( "A" ) letter in column 7 indicates that the triaxial deformation is (is not) included into the calculations. Column 
8 shows the pairing window used in the calculations with constant G and 8 pairing; no pairing window is used 
in the calculations with the Gogny force. Note that it is not always possible to extract these details from the 
original publications or references quoted therein. In these cases, the relevant box of the table is empty. 



Author [reference] 


Pairing Model 


Fitting region 


PNP 


Actinide 


SHE 


A/T 


-Bcutoff 


1 


2 


3 


4 


5 


6 


7 


8 


Macroscopic+microscopic method 


Moller 2009 [7] 


BCS(G)[48J 




Yes 


Yes 


Yes 


T 




Dobrowolski 2007 [6] 


BCS(G) 




No 


Yes 


2 


T 




Kowal 2010 [9] 


BCS(G) 


Z > 84 [51] 


No 


Yes 


Yes 


T 




Extended Thomas-Fermi plus Strutinsky integral 


Dutta 2000 [101 


BCS(<5) 




No 


5 


5 


T 




Skyrme density functional theory 


Bonneau 2004 [12] 


BCS(G)/BCS(<S) 


254 No/A ~ 178 


No 


Yes 


No 


A/T" 


6 MeV 


Biirvenich 2004 [52] 


BCS(5) 


across nuclear 
chart [53] 


No 


Yes 


Yes 


A 




Samyn 2005 [54] 


HFB (5) 


Z £ (92,98) c 


Yes 


Yes 


Yes 


T 


different E cutoa d 


Staszczak 2006 [13] 


BCS(G) 


282 Fm and Ref.[55j 


No 


Yes 


Yes 


T 


lowest Z(N) states' 3 


Staszczak 2007 [56] 


BCS(G)/BCS(<S) 


252 Fm 


No 


Yes 


Yes 


T 


lowest Z(N) states 6 


Gogny density functional theory 


Warda 2002 [3] 


HFB(Gogny) 


No 


No 


5 


No 


T 


No 


Delaroche 2006 [14] 


HFB (Gogny) 


No 


No 


Yes 


No 


T 


No 


Covariant density functional theory 


Bender 1998 [11] 


BCS(<5) 


across nuclear 
chart [53] 


No 


No 


3 f 


T 


[53/ 


Biirvenich 2004 [52] 


BCS(<5) 


across nuclear 
chart [53] 


No 


Yes 


Yes 


A 


[53/ 


Karatzikos 2010 [571 


RHB(Gogny) 9 


No 


No 


Yes 


Yes 


A 


No 


Abusara 2010 [1] 


BCS(G) 


Z € (90, 100) 
N- Z G (42,66) 


No 


Yes 


No 


T 


120 MeV 



"The calculations with allowance for triaxial deformation have been performed only for four nuclei 

b The pairing-active space Q, q is chosen to include approximately one additional oscillator shell of states above the Fermi 
level. 

Tairing force fitted to absolute masses. As a result, it is considerably stronger than the one fitted to even-odd mass 
differences [54j . 

d The single-particle states in the pairing window ep ± -E cu toff are included. -E cu tofT = 17 MeV for the BSk6, BSk7 and 
BSk8 Skyrme forces, E CVLtoB = 15 MeV for BSk2 force, E cutofi = 16.5 MeV for BSk9 force, E cutoS = 5 MeV for SLy6 force. 
e The pairing-active space consisted of the lowest Z(N) proton (neutron) single-particle states. 
^The calculations with allowance for triaxial deformation have been performed only for a single nucleus 
9 RHB=R,clativistic Hartrec Bogoliubov approach 
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FIG. 4: Potential energy surfaces of selected nuclei. The energy difference between two neighboring equipotential lines is equal 
to 0.5 MeV. The saddles along the 'Ax', 'Tr-A', and 'Tr-B' fission pathes are shown by solid circles, triangles, and squares, 
respectively. The solid diamonds show the outer fission barrier saddles. The saddles are defined via the immersion method 0, 
while the fission pathes as minimum energy pathes [HI, [64| which represent the most probable pathway connecting two minima 
via a given saddle. Note that contrary to other nuclei there are two possible fission pathes for outer fission barrier in the 300 116 
nucleus. 
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P 2 - deformation 



FIG. 5: (Color online) Deformation energy curves for 
the Z = 120, TV = 172 nucleus obtained with the NL3* 
parametrization. The black solid, red dashed and blue dot- 
dashed lines display the deformation energy curves for the 
axially symmetric, triaxial and axial octupole deformed solu- 
tions. We show the deformation energy curves for the last two 
solutions only in the range of /?2 values where it is lower in 
energy than the deformation energy curve of the axially sym- 
metric solution. Note that the 'Tr-B' fission path via saddle 
at /?2 = 0.47, 7 = 26° (see text for details) is not shown since 
this saddle is lower than the one of the 'Ax' fission path only 
by 0.2 MeV (see TableEJ ; the results of Ref. suggests that 
spontaneous fission exploiting this path is less likely than the 
one along the 'Ax' fission path. 
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N=172 N=174 N=176 N = 178 




P 2 deformation 

FIG. 6: (Color online) Deformation energy curves for the Z = 112, 114 and 116 nuclei obtained with the NL3* and DD-PC1 
parameterizations of the RMF Lagrangian. Solid lines correspond to axial solutions with reflection symmetry (A) , dashed lines 
to triaxial solutions with reflection symmetry (T), and dotted lines to octupole deformed solutions with axial symmetry (O). 
Note that the T and O solutions are shown only in the deformation range in which they are lower in energy than axial solution. 



N= 172 



N=174 N=176 N=178 




P 2 deformation 



FIG. 7: (Color online) The same as Fig. [6]but for the NL3* and DD-ME2 parametrizations. 



18 




P 2 cos(y+30) 

FIG. 8: The same as in Fig. [4] but for the nucleus 294 116. The results are shown for the parametrizations NL3*, DD-PC1 and 
DD-ME2. The energy difference between two neighboring equipotential lines is equal to 0.2 MeV. 




03 I I ■ I ■ I ■ I ■ I ■ I ■ I I 

164 168 172 176 180 184 188 
Neutron number N 

FIG. 9: (Color online) (panel (a)) The deformations of the 
ground states of even-even Z — 112 — 120 nuclei as a function 
of neutron number, (panel (b)) The energies of axial (large 
solid symbols, thick lines) and triaxial (either 'Ax-Tr' or 'Tr- 
A'; shown by small open symbols and thin lines) saddles with 
respect to the spherical/weakly (normal)-deformed minima in 
these nuclei as a function of neutron number. Note that the 
same type of symbols is used for axial and triaxial saddles in 
a given isotope chain. 
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TABLE V: The heights of axially symmetric ('Ax') and tri- 
axial ('Ax- TV', 'Tr-A' and 'Tr-B') saddle points [in MeV] 
with respect to spherical/weakly(normal)-deformed minima 
and their deformations. Note that only quadrupole deforma- 
tion is given for axial saddle. Columns 4 and 5 show the values 
for the saddles of either 'Ax-Tr' or 'Tr-A' fission pathes. The 
asterisk to the fission barrier height in column 4 indicates that 
the values for the 'Ax-Tr' path are displayed; the absence of 
the asterisk implies that the values for the 'Tr-A' path are 
shown. Note that triaxial saddles are shown only in the cases 
when their heights are lower than those of axial saddle. 



Nucleus 


Bf 


Pa 


Bf 


(P2,7) 


B jr- B 


Gfc.7) 


1 


2 


3 


4 


5 


6 


7 


Z = 120 nuclei 


308 120 
306 120 
304 120 
302 120 
300 120 
298 120 
296 120 
294 120 
292 120 


5.11 
5.00 
5.46 
5.37 
5.35 
5.35 
5.67 
5.27 
5.82 


0.20 
0.20 
0.21 
0.19 
0.32 
0.33 
0.34 
0.34 
0.36 


3.04 
3.46 
4.51 
5.00 
4.93* 
5.02* 
5.57* 


(0.34,26°) 
(0.33,26°) 
(0.33,26°) 
(0.33,25°) 
(0.31,11°) 
(0.33,12°) 
(0.35,10°) 


2.23 
2.02 
3.30 
3.93 
4.08 
4.27 
4.92 
4.73 
5.61 


(0.48,22°) 
(0.50,22°) 
(0.51,21°) 
(0.50,24°) 
(0.50,26°) 
(0.49,28°) 
(0.50,27°) 
(0.50,27°) 
(0.50,27°) 


Z = 118 nuclei 


304 118 
302 118 
300 118 
298 118 
296 118 
294 118 
292 118 
290 118 

288 n8 


4.97 
5.07 
5.18 
5.07 
5.00 
4.95 
4.75 
4.15 
3.28 


0.21 
0.23 
0.28 
0.32 
0.32 
0.32 
0.34 
0.36 
0.36 


3.92 
4.57 
4.94* 
4.72* 
4.73* 
4.82* 


(0.33,26°) 
(0.33,26°) 
(0.27,15°) 
(0.32,12°) 
(0.33,10°) 
(0.33,10°) 


3.03 
3.41 
3.37 
3.67 
3.87 
4.12 
4.35 
4.08 


(0.51,26°) 
(0.51,25°) 
(0.51,26°) 
(0.50,27°) 
(0.50,26°) 
(0.50,25°) 
(0.50,23°) 
(0.49,24°) 


Z — 116 nuclei 


300 116 
298 116 
296 116 
294 116 
292 116 
290 116 
288 116 
286 116 
284 116 


4.83 
4.78 
4.53 
4.37 
3.98 
3.39 
3.55 
3.43 
3.40 


0.26 
0.28 
0.32 
0.32 
0.33 
0.34 
0.34 
0.36 
0.37 


4.48 
4.57 
4.32* 
4.22* 
3.93* 


(0.35,25°) 
(0.30,15°) 
(0.32,11°) 
(0.33,7°) 
(0.34,7°) 


3.13 
3.29 
3.35 
3.40 
3.39 
3.47 


(0.52,22°) 
(0.51,23°) 
(0.51,24°) 
(0.50,25°) 
(0.48,26°) 
(0.48,27°) 


Z = 114 nuclei 


296 114 
294 114 
292 114 
290 114 
288 114 

286 U4 

284 114 
282 114 

280 1U 


4.65 
4.19 
3.81 
3.40 
3.52 
3.32 
3.07 
3.08 
3.20 


0.28 
0.29 
0.32 
0.34 
0.32 
0.33 
0.38 
0.42 
0.42 


4.43* 
4.14* 


(0.31,15°) 
(0.31,10°) 


3.42 
3.08 
3.18 
3.38 


(0.49,24°) 
(0.49,25°) 
(0.48,25°) 
(0.48,26°) 
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TABLE VI: The same as Table M but for Z = 112 nuclei. 



Nucleus B 



Bf 



Cfti,7) 



7 



112 nuclei 



292 112 


3.94 


0.30 


290 112 


3.48 


0.31 


288 112 


3.07 


0.31 


286 112 


2.90 


0.32 


284 112 


2.65 


0.32 


282 112 


2.35 


0.41 


280 112 


2.70 


0.43 


278 112 


3.06 


0.44 


276 112 


3.34 


0.44 



3.64 
3.44 



(0.48,27°) 
(0.47,26°) 



TABLE VII: The same as in Table fVl but only for nuclei where 
experimental estimates of inner fission barrier heights exist. 
The results of calculations with the DD-PC1 and DD-ME2 
parametrizations are presented. 



Nucleus 


Bf 


P2 


Bf 


C&, 7 ) 


B T r B 


G&.7) 


DD-PC1 parametrization 


294 116 


3.64 


0.31 


3.47* 


(0.32,10°) 


2.58 


(0.48,27°) 


292 U6 


3.63 


0.32 


3.55* 


(0.32,5°) 


2.83 


(0.49,26°) 


290 114 


3.42 


0.34 










288 U4 


3.38 


0.32 










286 U2 


3.02 


0.31 










284 U2 


2.75 


0.31 










DD-ME2 parametrization 


294 n6 


5.36 


0.32 


5.20* 


(0.32,8°) 


4.45 


(0.48,25°) 


292 U6 


4.77 


0.32 


4.68* 


(0.32,5°) 


4.24 


(0.48,27°) 


290 114 


4.33 


0.34 










288 U4 


4.18 


0.32 










286 U2 


3.62 


0.31 










284 112 


3.12 


0.31 
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FIG. 10: (Color online) The heights of outer fission bar- 
riers of even-even Z — 112 — 120 nuclei relative to 
spherical/weakly(normal)-deformed minima as a function of 
neutron number. The results of calculations of type A (ax- 
ial), T (triaxial) and O (octupole) obtained with the NL3* 
parametrization are presented. 
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1.25 



deformation 

FIG. 11: (Color online) The same as in Fig. [5] but for the 
nuclei 298 120 and 306 1 20. Note that for simplicity only the 
axial solution (A) is shown at ft? < 0.5. 




0.25 0.5 

P 2 cos(y+30) 



FIG. 12: (Color online) The same as in Fig. [4] but for the nu- 
cleus 302 120. In order to show the potential energy surface be- 
tween the axial /3 ~ 0.2 and triaxial (p% ~ 0.40, 7 ~ 25°) hills 
in greater details the energy difference between two neighbor- 
ing equipotential lines is set to 0.1 MeV and only one fourth 
of deformation plane of Fig. [4] is shown. 
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TABLE VIII: A summary of the nuclei in which the superde- 
formed minimum is the lowest in energy. The energies of outer 
fission barrier saddle points with respect to this minimum, 
as obtained in axial (Bf_ n ), triaxial (Bj_ II ) and octupole 
(Bf_jj) deformed calculations, are shown in columns 3, 4, 
and 5, respectively. For each nucleus, the outer saddle point 
with lowest energy is shown in bold. The excitation energies 
E exc of spherical/ weakly (normal) -deformed minima with re- 
spect to the superdeformed minima are shown in column 2. 
The graphical explanation of these quantities is also given in 

Fig.rru 



Nucleus 




Bf-n 


Bj-ii 


Bf-H 


1 


2 


3 


4 


5 


The Z = 120 nuclei 


im l2Q 


4.84 


4.52 


2.75 


1.62 


3UB 120 


4.26 


4.36 


3.00 


1.82 


a,J4 120 


3.48 


4.49 


3.57 


2.07 


3U *120 


2.89 


5.30 


3.65 


2.55 


auu 120 


2.60 


6.18 


3.82 


2.61 


!«B 120 


1.99 


5.71 


3.38 


2.25 




1.72 


4.87 


2.94 


2.52 


«* 120 


0.62 


4.63 


2.24 


2.71 


The Z = 118 nuclei 


3U4 llg 


2.68 


2.41 


2.14 


1.31 


aua 118 


2.92 


3.58 


2.42 


1.42 


a.)u 118 


1.78 


5.16 


3.40 


2.20 


!1UB 118 


1.70 


5.34 


3.54 


2.34 


m 118 


1.48 


5.17 


3.15 


2.26 


!«4 118 


1.15 


4.99 


2.82 


2.35 


^118 


0.69 


4.58 


2.28 


2.59 


ayu 118 


0.40 


4.03 


1.61 


2.60 


The Z = 116 nuclei 


aJU 116 


0.92 


2.76 


2.60 


1.75 


m 116 


0.83 


3.98 


2.95 


1.88 



